The differential information-geometry of quantum phase transitions 
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The manifold of coupling constants parametrizing a quantum Hamiltonian is equipped with a 
natural Riemannian metric with an operational distinguishability content. We argue that the sin- 
gularities of this metric are in correspondence with the quantum phase transitions featured by the 
corresponding system. This approach provides a universal conceptual framework to study quantum 
critical phenomena which is differential-geometric and information-theoretic at the same time. 
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Introduction- Suppose you are given with a set of 
quantum states associated to a family of Hamiltonians 
smoothly depending on a set of parameters, e.g., coupling 
constants. This parameter manifold - that can include 
temperature in the case the considered states are thermal 
ones - is partitioned in regions characterized by the fact 
that inside them one can "adiabatically" move from one 
point to the other and no singularities in the expectation 
values of any observables are encountered. The bound- 
aries between these regular regions are in turn associated 
to the non-analytic behaviour of some observable and are 
referred to as critical points; crossing one of these points 
results in a phase transition (PT). States lying in different 
regions generally have some strong structural difference 
and are, in principle, easily distinguishable once somehow 
a preferred observable is chosen. 

The standard machinery, i.e., the so-called Landau- 
Ginzburg paradigm, to deal with this phenomenon is 
based on the notions of symmetry breaking, order pa- 
rameter, and correlation length On the other hand, 
some system fails to fall in this conceptual framework. 
This can be due to the difficulty of identifying the proper 
order parameter for systems whose symmetry breaking 
pattern is unknown or to the very absence of a local or- 
der parameter, e.g., quantum phase transitions (QPTs) 
involving different kinds of topological order Even 
another standard characterization of QPTs, i.e., singu- 
larities in the ground state (GS) energy as a function of 
the coupling constant, misses to capture the boundaries 
between phases for some QPTs, e.g., those with matrix- 
product states 

In the last few years ideas and tools borrowed from 
quantum information science Q have been used to study 
quantum, i.e., zero temperature, phase transistions [5(; 
in particular the role of quantum entanglement in QPTs 
has been extensively investigated [f|. More recently an 
approach to QPTs based on the concept of quantum fi- 
delity has been put forward 0] and applied to systems of 
quasi-free fermions [1, E| , to the so-called matrix-product 
states and extended to finite-temperature HI]. ^ n 
the fidelity approach, QPTs are identified by studying 



the behavior of the amplitude of the overlap, i.e., scalar 
product, between two ground states corresponding to two 
slightly different set of parameters; at QPTs a drop of 
the fidelity with scaling behaviour is observed and quan- 
titative information about critical exponents can be ex- 
tracted [I, [l(| • The fidelity approach is not based on the 
identification of an order parameter - and therefore does 
not require a knowledge of symmetry breaking patterns 
- or more in general on the analysis of any distinguished 
observable, e.g., Hamiltonian, but it is a purely metrical 
one. All the possible observables are in a sense considered 
at once. 

In this paper we shall unveil the universal differential- 
geometric structure underlying these observations. We 
shall show how QPTs can be associated to the singu- 
larities of a Riemannian metric tensor inherited by the 
parameter space from the natural Riemannian structure 
of the projective space of quantum states. This structure 
has an interpretation in terms of information-geometry 
[H, EH providing the differential-geometric approach of 
this paper with an information-theoretic content. 

Information- geometry and QPTs- Let us consider a 
smooth family H(\), A 6 A4(=the parameter manifold), 
of quantum Hamiltonians in the Hilbert-space TC of the 
system. If |^ (A)) € H denotes the (unique for sim- 
plicity) ground-state of H(X), one has defined the map 
^o : M — > H/X — ► |^o(A)) associating to each set of 
parameters the ground-state of the corresponding quan- 
tum Hamiltonian. This map can be seen also as a map 
between M. and the projective space P7i(=manifold of 
"rays" of Ti). This space is a metric space being equipped 
with the so-called Fubini-Study distance dFs{i } , < t>) : = 
cos^ 1 F(ip, <fi), where 



J^M) := KiM>l 



(1) 



and 1 1 ip\ j = || 0| I = 1. In Ref. [12( Wootters showed that 
this metric has a deep operational meaning: it quanti- 
fies the maximum amount of statistical distinguishability 
between the pure quantum states \ip) and \<f>). More pre- 
cisely, dFsi^i 4>) is the maximum over all possible projec- 
tive measurements of the Fisher-Rao statistical distance 



2 



between the probability distributions obtained from 
and \<j)) Moreover, this result extends to mixed 

states as well by replacing the pure-state fidelity (Q]) with 
the Uhlmann fidelity 15| and the projective measure- 
ments with generalized ones [l3j |. 

These results are non-trivial and allow to identifty in 
a precise manner the Hilbert space geometry with a ge- 
ometry in the information space: the bigger the Hilbert 
(or projective) space distance between \ip) and \<f>) the 
higher the degree of statistical distinguishability of these 
two states. From this perspective it is clear that a single 
real number, i.e., the distance, virtually encodes informa- 
tion about infinitely many observables, e.g., order param- 
eters, one may think to measure. This remark basically 
contains the main intuition at the basis of the metric 
approach to QPTs advocated in this paper: at the tran- 
sition points, a small difference between the control pa- 
rameters results in a greatly enhanced distinguishability 
of the corresponding GSs, which should be quantitatively 
revelead by the behavior of their distance. 

For the purposes of this paper it is crucial to note 
that the projective manifold PH, besides the structure of 
metric space, has a well-known structure of Riemannian 
manifold, i.e., it is equipped with a metric tensor. Here, 
for the sake of self-consistency, we briefly recall how this 
Riemannian metric is obtained starting from the Hilbert 
space structure of H. PH can be seen as the base man- 
ifold of a (principal) fiber bundle with total space given 
by the unit ball S of H, i.e., S := {\ip) E H / \\ip\\ = 1}, 
and projection tt: S — > PH / \ip) — » {e ie \ip) / 6 e [0, 2tt)}. 
The tangent space to each point of S is isomorphic 
to a subspace of H and has therefore defined over it 
the Hermitean bilinear form g\^\(u, v) := (u, v) (u and 
v are tangent vectors, i.e., elements of H). This defines 
a (complex) metric tensor field g over S. To project g 
down to PH one has to introduce the notion of hori- 
zontal subspace for each tangent space of S or equiva- 
lently that of parallel transport and the associated one 
of connection. In this case the Hilbert space struc- 
ture of the tangent spaces provides a natural solution 
to this task: the horizontal subspace is simply the set 
of vectors \u) which are orthogonal to the fiber over \i/j), 
i.e., {u,ip) = 0. It follows that the complex metric over 
PH is given by g n{W) (u,v) = (ujl - \i>) (ip\)v) , called 
the quantum geometric tensor [la ]. The real (imagi- 
nary) part of this quantity defines a Riemannian metric 
tensor (symplectic form) on PH. Another, elementary 
way of getting the form of the Riemannian metric over 
PH is by means of Eq. ([T]). For T very close to the 
unity, one can write d 2 FS {i\),i\) + Sip) — 2(1 — J 7 ). Since 
^(tp^ + Stp) ~ \l + (ip,5ip) + (l/2){iP,S 2 ip)\ 2 , using this 
expression and the normalization of one finds 

ds 2 -. = 4 s (^+ty) = w,^)-K^^)r 

= (^(l-hMviW) • (2) 



the parameter manifold M. induced, i.e., "pulled-back" 
by the ground state mapping Wo introduced above. By 
writing J|tf (A)) = E^MdX", with 9 M := 0/flA", 
fi = 1, . . . , dim M., and using Eq. ((2|), one imediately 
obtains ds 2 = g^ dX^dX", where 



(3) 



Now we provide a simple perturbative argument on why 
one should expect a singular behavior of the metric ten- 
sor at QPTs [l8|. By using the first order perturba- 
tive expansion |* (A + 8X)) ~ |*o(A)> + E^oC^o - 
S„)- 1 |*„(A))(*„(A)|^|*o(A)), where SH := H{X + 
SX) — H(X), one obtains for the entries of the metric ten- 
sor ([3]) the following expression 



9^ 
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(Eo (A)|0 M ff|tt n (A))(tt n (A)|0,,ff|tt o (A)) 
[E n (X) - E (X)} 2 



What we would like to do now is to see the metric in 



(4) 

An analogous expression, with the real part replaced by 
the imaginary one, gives the antisymmetric tensor which 
describes the curvature two-form whose holonomy is the 
Berry phase [17[ ■ Continuous QPTs are known to occur 
when, for some specific values of the parameters and in 
the thermodynamical limit, the energy gap above the GS 
closes. This amounts to a vanishing denominator in Eq. 
(j4| that may break down the analyticity of the metric 
tensor entries. 

To get further insight about the physical origin of these 
singularities we notice that the metric tensor Q can be 
cast in an interesting covariance matrix form [la] . In 
the generic case, by moving from H(X) to H(X + 5X) 
no level-crossings occur. In this case the unitary opera- 
tor 0(X,5X) := E„l*n(A + 5A)(# n (A)| "adiabatically" 
maps the eigenvectors at A onto those at A + <5A. Then by 
introducing the observables := i(d fl O)0^ the metric 
tensor ^ takes the form g^ = (1/2)({A > M , X v }) where 
Xp := X^ — (Xp). Moreover, the line element ds 2 can 
be seen as the variance of the observable X := i(dO)0\ 
i.e., ds 2 = (X 2 ). The operator X is the generator of the 
map transforming eigenstates corresponding to different 
values of the parameter into each other. The smaller the 
difference between these eigenstates for a given parame- 
ter variation, the smaller the variance of X. Intuitively, 
at the QPT one expects to have the maximal possible 
difference between |5 , (A)) and |^o(A + SX)), i.e., many 
"unperturbed" eigenstates |^ n (A)) are needed to build 
up the "new" GS; accordingly the variance of X can get 
very large, possibly divergent. In a sense ds 2 can be seen 
as a sort of susceptibility of the "order parameter" X. 

Quasi-Free fermionic systems.- In order to show ex- 
plicitly how the singularities, i.e., divergencies of g^ u 
arise, we will discuss the case of the XY model in a de- 
tailed fashion; before doing that we would like to make 
some general considerations about the systems of quasi- 
free fermions on the basis of the results presented in Ref. 
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[8j. Systems of quasi- free fermions are defined by the 
following quadratic Hamiltonian 



l ^ L 



(5) 



'■j=i 



where: the Cj's (cj's) are the annihilation (creation) 
operators of L fermionic modes, A, B G M_(R) are 
L x L real matrices, symmetric and anti-symmetric re- 
spectively, i.e., A T = A, B T = — £?. In Ref. @| 
it has been shown that the set of GSs of Eq. ((Sj) 
is parametrized by orthogonal L x L real matrices T 
giving the unitary part of the polar decomposition of 
the matrix Z := A — B. One can then prove that 
T{Z,Z') := |(*z|* Z /)| = V|det[(r + T')/2]| [8]. With 
no loss of generality we can assume det(T) = 1 which 
identifies the GS manifold of the quasi-free systems ((Sj) 
with SO(L,R). Since /(_"') := F(Z, Z') has a maximum 
equal to one at Z' = Z one has S 2 f(Z')\ z = <5 2 In /(_"') |_; 
from this, the expansion for Z' — > Z of the above for- 
mula for T (Eq. (8) in Ref. Q) and by defining 
K := InT € so(i, R), one finds an explicit form for the 
metric: ds 2 ~ 2(1 - T) = (l/8)Tr(dAT) 2 . From this equa- 
tion, if K = K(X), with A e M, one obtains the fol- 
lowing expression for the metric tensor induced over M., 
i-e., g^v = (l/8)Tr(d IJr Kd v K). For translationally invari- 
ant Hamiltonians (JSJ) the anti-symmetric matrix if can 
be always cast in the canonical form K = i®k@kO~k where 
k is a momentum label. Therefore in this important case 
one has fl „„ = (l/4)£ fc (_0 fc /0A")(_0fc/_A''). 

We see here that the connection established in Refs. 
0) between QPTs, e.g., due to the vanishing of a 
quasi-particle energy, and a singularity in the second or- 
der expansion of T can be directly read as a connection 
between QPTs in quasi-free systems and singularities in 
the metric tensor g^ v . 

The nature of this connection will be now exempli- 
fied by considering the QPTs of the periodic antiferro- 
magnetic XY spin chain in a transverse magnetic field. 
By writing the spin operator in terms of Pauli matri- 
ces, i.e., S = it/2, the Hamiltonian for an odd num- 



ber of spins L = 2M + 1 reads H 



Ml 



"(I 



7)ojfff+i/4 - (1 - 7)f"o-"+i/4 + ha]/2], where 7 is the 
anisotropy parameter in the x-y plane and h is the mag- 
netic field. This Hamiltonian can be cast in the form 
([5]) by the Jordan- Wigner transformation. The criti- 
cal points of this model are given by the lines h = ±1 
and by the segment \h\ < 1,7 = 0. The single parti- 



cle energies are = J e| + 7 2 sin {2irk/L), where ek = 
cos(2irk/L) — h and k = —M, . . . , M. For this model the 
Ok 's defined above have the form 9k = cos _1 (efe/Afc) and 

<W = (±/4)Ek=i(d0k/d\n(d9k/d\n, where A 1,2 = 
fc,7. One finds (d6 k /dh) 2 = 7 2 sin 2 x k /K\, (d6 k /dj) 2 = 
sin 2 Xk(cosxk — h) 2 / 'A*, and (d8k/dh)(d8k/dj) — 
7 sin 2 Xk (cos Xfe — h)/K\, with 2^ = 2nk/L. 




FIG. 1: Induced curvature _R scaled by the system size L for 
the parameter space of the XY model. 



In the thermodynamic limit (TDL), the explicit cal- 
culation of g^ u can be performed analytically. Indeed, 
except at critical points, for large L one can replace 
the discrete variable Xk with a continuous variable x 
and substitute the sum with an integral, i.e., 53fe=i ~* 
[L/(2-k)\ J dx. At critical points this is not generally 
feasible due to singularities in some of the terms in the 
sums. Outside critical points, the resulting integrals, al- 
beit non-trivial, yield simple analytical formulas, which 
differ depending on whether \h\ < 1 or \h\ > 1. 

For \h\ < 1 in the TDL one finds a diagonal metric 
tensor 



9 



16| 7 



■diag 



1 



1 



1_^'(1 + | 7 |)5 



(6) 



Closed analytic formulas in the TDL can be obtained 
also for I h \ > 1, although in a less compact form, which 
we omit here for brevity. We only note that for \h\ > 1 
also the off-diagonal elements of the metric tensor are 
non-zero. Having the induced metric tensor it is also 
possible to investigate the induced curvature of the pa- 
rameter manifold. We therefore compute the scalar cur- 
vature R, which is the trace of the Ricci curvature tensor 
[H. We find R(\h\ < 1) = -(16/L)( 1 + J7D/I7I and 
R(\h\ > 1) = {16/L)(\h\ + V/i 2 +7 2 -l)/\A 2 + 7 2 -l- 
Note that the curvature diverges on the segment \h\ < 
1,7 = and is discontinuous on the lines h = ±1. In- 
deed, lhn|/ 1 ,|_ ) .i+ R = —limift.i_i.i-ii. The behaviour of 
the curvature R is shown in Fig. [T] 

Mixed states and classical transitions - In this section 
we would like to make some extentions of the idea de- 
veloped in this paper to finite temperature. This will 
allow us to establish a connection between the present 
approach and the one for classical PTs developed in 
20l 121 1 . This latter formalism is in fact obtained in the 
special case of commuting density matrices which effec- 
tively turns the quantum problem into a classical one. 

The fidelity approach to QPTs can be extended to 
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by using the 

r 1/2 l/2il/ 2 

TrK P0P1 In- 



finite-temperature, i.e., to mixed-states, 
Uhlmann fidelity [H|: T{p$,pi) 
When po and p± are commuting operators the fidelity 
takes the form T(po,p\) = Y~ \/p®vL where the 
are the eigenvalues of the p Q 's [23[. In particular, when 
p a = Z- 1 exp(- /3 a H),Z a := Trexp(-/3 a H), (a = 0,1) 
one immediately finds that the fidelity has a simple ex- 
pression in terms of partition functioms: T = Z((3 /2 + 



(3 1 /2)(Z(/3 )Z(f3 1 ))- 1 ^ 
(3,(3i — (3 + 5(5 one obtains 



T{(3, f3 + 6/3)~ exp 



111 ]. By expanding for /3 = 



J/? 2 
'8/3' 



(7) 



where cy (/?) denotes the specific heat [l| . This relation is 
remarkable in that it connects the distinguishability de- 
gree of two neighboring thermal quantum states directly 
to the macroscopic thermodynamical quantity cy. The 
line element of the parameter space, i.e., the (3 axis, is 
then given by ds 2 ~ c v (f3)P~ 2 dp 2 = {{H 2 ) p - {HfydJ 2 . 
A closely related formula has been obtained in 2Q, 21 1- 
Since PTs are associated to anomalies, e.g., divergences, 
in the behavior of cy (/?) , we see that also in this "classi- 
cal" case the metric ds 2 induced on the parameter space 
contains signatures of the critical points. In this sense 
the information-geometrical approach to PTs seems able 
to put quantum and classical PTs under the same con- 
ceptual umbrella. 

Conclusions - In this paper we proposed a differential- 
geometric approach to study quantum phase transtions. 
The basic idea is that, since distance between quantum 
states quantitatively encodes their degree of distinguisha- 
bility, crossing a critical point separating regions with 
structurally different phases should result in some sort of 
singular behaviour of the metric. This intuition, based 
on early studies of quantum fidelity, can be made rigor- 
ous in some simple yet important cases, e.g., quasi- free 
fermion systems. The manifold of coupling constants pa- 
rameterizing the system's Hamiltonian can be equipped 
with a (pseudo) Riemannian tensor g whose singularities 
correspond to the critical regions. For the case of the 
XY chain we explicitely computed the components of g 
in the thermodynamic limit, showing that they are diver- 
gent, with universal exponents, at the critical lines. We 
also computed the scalar curvature of g and analyzed its 
relation with criticality. The geometrical approach advo- 
cated in this paper does not depend on the knowledge of 
any order parameter or on the analysis of a distinguished 
observable, it is universal and information-theoretic in 
nature. The study of the physical meaning of the geo- 
metric invariants one can build starting from g (e.g., the 
curvature), their finite-size as well as scaling behaviour, 
and their relations with the nature of the quantum phase 
transition are important questions to be addressed in fu- 
ture research. 
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